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Abstract 



We regard drift-diffusion equations for semiconductor devices in Lebesgue 
spaces. To that end we reformulate the (generalized) van Roosbroeck sys- 
tem as an evolution equation for the potentials to the driving forces of the 
currents of electrons and holes. This evolution equation falls into a class of 
quasi-linear parabolic systems which allow unique, local in time solution in 
certain Lebesgue spaces. In particular, it turns out that the divergence of 
the electron and hole current is an integrable function. Hence, Gauss' theo- 
rem applies, and gives the foundation for space discretization of the equations 
by means of finite volume schemes. Moreover, the strong differentiability of 
the electron and hole density in time is constitutive for the implicit time dis- 
cretization scheme. Finite volume discretization of space, and implicit time 
discretization are accepted custom in engineering and scientific computing. 

- This investigation puts special emphasis on non-smooth spatial domains, 
mixed boundary conditions, and heterogeneous material compositions, as re- 
quired in electronic device simulation. 
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1 Introduction 

In 1950 van Roosbroeck [48] established a system of partial differential equations 
describing the motion of electrons and holes in a semiconductor device due to drift 
and diffusion within a self-consistent electrical field. In 1964 Gummel [28] published 
the first report on the numerical solution of these drift-diffusion equations for an op- 
erating semiconductor device. From that time on van Roosbroeck's system has been 
the backbone of many a model in semiconductor device simulation. The first papers 
devoted to the mathematical analysis of van Roosbroeck's system appeared in the 
early seventies of the previous century [38, 39]; for a historical synopsis and further 
references see [11]. In 1986 Gajewski and Groger proved the global existence and 
uniqueness of weak solutions under realistic physical and geometrical conditions 
[13]. The key for proving these results and also for establishing stable numerical 
solving procedures is the existence of a Lyapunov function for the van Roosbroeck 
system. This solution theory entails restricting conditions on the models for the 
recombination of electron-hole pairs, see [11, 2.2.3], [14, Ch. 5], [15, Ch. 6], [18], and 
[19]. In this paper we relax the condition on the reaction terms in the equations 
considerably, up to the point that some external control to the generation or anni- 
hilation of electrons or holes can be applied individually. In particular, this aims at 
radiative recombination of electron-hole pairs in semiconductor lasers, and at the 
generation of electron-hole pairs in optoelectronic detectors. Notwithstanding this 
generalization, we continue to use the name van Roosbroeck system for the model 
equations. 

Van Roosbroeck's system consists of current-continuity equations — one for elec- 
trons, another one for holes — which are coupled to a Poisson equation for the 
electrostatic potential, and comprise generative terms, first of all recombination of 
electron-hole pairs. The current-continuity equations can be viewed as quasi-linear 
parabolic equations. However, the natural formulation of balance laws is in integral 
form 



Here U2 and u\ is the density of electrons and holes, respectively, jk is the corre- 
sponding flux, and is a reaction term, uj is any (suitable) sub-domain of the 
whole domain under consideration, v the outer unit normal to the boundary du of 
uj and cr^ the arc measure on du. In the weak formulation of the balance law the 
boundary integral of the normal component of the current is expressed as the volume 
integral of the divergence of the corresponding current. Very little is known about 
the question whether the weak solutions also satisfy the original balance law equa- 
tions (1.1). Obviously, this depends on the applicability of Gauss' theorem. So, the 
problem is about the divergence of the currents in weak solutions being functions - 




(1.1) 
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not only distributions. In particular, this comes to bear in the numerical treatment 
of van Roosbroeck's system. The choice for space discretization of drift-diffusion 
equations is the finite volume method, see [17], which rests on the original balance 
law formulation (1.1) of the equations. 

In this paper we solve this problem for the spatially two-dimensional van Roosbroeck 
system by showing that it admits a classical solution in a suitably chosen Lebesgue 
space — at least locally in time. Aiming at the inclusion of rather general recom- 
bination and generation processes for electron-hole pairs we cannot expect global 
existence anymore, and we cannot rely on a Lyapunov function. Instead we apply 
local methods for quasi-linear evolution equations. To that end, we rewrite van 
Roosbroeck's system as an evolution equation for the electrochemical potentials of 
electrons and holes, and apply a recently obtained result on quasi-linear parabolic 
equations in Lebesgue spaces, see [31]. This yields a classical solution of van Roos- 
broeck system locally in time with currents the divergence of which is Lebesgue 
integrable to some exponent greater than one. The strong differentiability of the 
electron and hole density in time is constitutive for the implicit time discretization 
scheme which is accepted custom in engineering and scientific computing, see for 
instance [11]. 

Please note that in device simulation one is always confronted with contacted devices 
of heterogeneous material composition. That leads to mixed boundary conditions 
and jumping material coefficients in the model equations. Hence, standard theorems 
on existence, uniqueness and regularity do not apply. 



2 Van Roosbroeck's system 

Basic variables 

In the following we investigate van Roosbroeck's model for a semiconductor device 
which describes the flow of electrons and holes in a self-consistent electrical field due 
to drift and diffusion. The physical quantities one is interested in are: the densities 
Mi and U2 of holes and electrons, the densities ji and ji of the hole and electron 
current, the electrostatic potential tp of the self-consistent electrical field, and the 
electrochemical potentials <f>\ and 02 of holes and electrons These unknowns have 
to satisfy Poisson's equation and the current-continuity equations for electrons and 
holes with some side conditions. The latter are given by the relations between the 
potentials and the densities. 
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Spatial domain 

We study only semiconductor devices which are quasi translational invariant in one 
space direction or angular symmetric. In that case van Roosbroeck's system in real 
space can be reduced to a similar set of equations in the plane. That means, we 
regard a cut through the device perpendicular to the direction of invariance. Let 
Q be the resulting two-dimensional (bounded) representative domain. Parts of the 
device may be insulating, for instance formed by an oxide. Then, electrons and holes 
can move only in a sub-domain Q of Q. This also covers the case of charges which 
are artificially immobilized on a sub-domain Furthermore, we mark out a part 

T of the boundary of Q where the device borders on an insulator. The remaining 
part of the boundary represents (possibly several) contacts of the device. We also 
mark out a part T of fi's boundary. In the case of a stand alone drift-diffusion 
model of the semiconductor device again T represents areas of the device bordering 
to an insulator, whereas the remaining part is the contact area. 

External control 

In real-world modeling of semiconductor devices van Roosbroeck's system often 
serves as a component in a compound model of the device. Then the superordinated 
system — for instance a circuit model — may exercise a control on van Roosbroeck's 
system. Apart of a superordinated circuit model, compound models comprising in 
addition to van Roosbroeck's system equations for the lattice temperature or the 
power of lasing modes play an important role in device simulation, see for instance 
[11, 2, 4, 3]. But the concept of external control also comes to bear in segmentation 
of the simulation domain, in particular in connection with multiscale modeling, see 
for instance [32, 33, 30]. 

If van Roosbroeck's equations serve as a component of a compound model, then 
system parameters, state equations, boundary conditions, et alii, possibly bear a 
different physical meaning than in the stand-alone model. 

We make assumptions about an external control from the initial time T up to a 
time T]_. 

2.1 Poisson equation 

The solution of the Poisson equation with mixed boundary conditions, 




<P = Vd(*) 






on r, 
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gives the electrostatic potential tp on ft subject to the electron and hole density u 2 
and u\. Strictly speaking, the densities u k , k = 1,2, are only defined on ft but, we 
extend them by zero to ft. 

The parameters in (2.1) have the following meaning: e is a bounded, measurable 
function on ft with values in the set of real, symmetric, 2x2, positive definite 
matrices and corresponds to the spatially varying dielectric permittivity on the space 
region occupied by the device. Moreover, we assume 

lk(^)||B(R 2 ) < and ( e ( x )0 ' £ > £ «II£||r2 f° r almost all x e ft and all £ G IR 2 

with two strictly positive constants e. and e*. Furthermore, is a non-negative 
function on T, representing the capacity of the part of the device surface bordering 
on an insulator. We assume that D is not empty or ep is positive on a subset of 
T with positive arc measure. In other words, the device has a Dirichlet contact or 
part of its surface has a positive capacity. y?g(t) an d V 9 f(^) are ^ ne v °hages applied 
at the contacts of the device, and d(t) represents a charge. In the case of a stand 
alone drift-diffusion model <pq, <pf, and d are constant in time, and d solely is the 
charge density of dopants in the semiconductor materials composing the device. In 
general, <pq, <pf, and d are function which are defined on the time interval [To,Xi] 
where a possible control acts on the device. 



2.2 Current— continuity equations 

The current -continuity equations for holes and electrons (k = 1,2, respectively) 

u' k - V • j k = r k (t, <p, 4> x , 4> 2 ) on ft (2.2) 

characterize the evolution of the electron and hole density under the action of the 
currents j k and the reactions r k subject to the mixed boundary conditions 

<t>k{t) = 4>T>,k{t) on D = interior(9ft \ T), ^ 

vj k = on T, 

from the initial conditions 

MT ) = $° fc - (2.4) 

Each r k , k — 1, 2 is a reaction term which models the generation and annihila- 
tion of electrons and holes. In particular, this term covers the recombination of 
electrons and holes in the semiconductor device. r\ and r 2 can be rather general 
functions of the particle and current densities, see §2.4. We require that the set 
D = interior (9ft \ T) is not empty. The boundary values 0d,i, 0d,2 in general de- 
pend on time. Moreover, the reactions r k may explicitly depend on time. This 
dependence on time, again, allows for a control of the system by some other part of 
a superordinated compound model. 
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2.3 Carrier and current densities 

Van Roosbroeck's system has to be complemented by a prescription relating the 
density of electrons and holes as well as the densities of the electron and hole current 
to the chemical potentials of these charge carriers. We assume 

u k (t,x) = p k (t,x)J 7 k (xk{t,x)), xett, k = 1,2, (2.5) 

where Xi and X2 are the chemical potentials 

Xk = 0k + {-l) k V + b k , k = 1,2, (2.6) 

and 02, 0i are the electrochemical potentials of electrons and holes, respectively. 
b k , Pk, k = 1,2 are positive, bounded functions on Q. They describe the electronic 
properties of the materials composing the device. b 2 and b\ are the band edge offsets 
for electrons and holes, and p 2 , pi are the corresponding effective band edge densities 
of states. If the equations under consideration form part of a compound model for 
the semiconductor device, then b k , p k , k = 1,2, may depend on time. For instance, 
the p k could be subject to an external control of the device temperature. Then they 
depend on time via the temperature. Mathematically, we assume the following. 

2.1 Assumption. For every t e [Tq,Ti] the functions p k (t) are essentially bounded 
on Q and admit positive lower bounds which are uniform in t G [T ,Ti]. The 
mappings 

[T , TJ 3 t h-> p k {t) EL 2 (n), k = 1,2 (2.7) 
are different iable on the interval ]T ,Ti[ with Holder continuous derivatives p' k . 

The functions T\ and T 2 represent the statistical distribution of the holes and elec- 
trons on the energy band. In general, Fermi-Dirac statistics applies, i.e. 

Ms) = 4= r tt^i ^ seR - ( 2 - g ) 

V7r J 1 + e* 

However, often Boltzmann statistics F^s) = e s is a good approximation. 

As for the kinetic relations specifying the current-continuity equations we assume 
that the electron and hole current is driven by the negative gradient of the electro- 
chemical potential of electrons and holes, respectively. More precisely, the current 
densities are given by 

j k {t,x) = -Q k {xk{t,x)) p k {x)V4> k {t,x) , 16S], k = l,2. (2.9) 

The mobilities p 2 and p\ for the electrons and holes, respectively, are measurable, 
bounded function on Q with values in the set of real, 2x2, positive definite matrices 
satisfying for almost all x G Q and all (el 2 

HMaOIki; 2 ) < At* and (jik(?)€) ■ £ > A*.||£||r2, k = 1,2, 
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with two strictly positive constants /i. and /i*. The mobilities are accounted for on 
the parts of the device where electrons and holes can move due to drift and diffusion. 

2.2 Remark. In semiconductor device modeling, usually, the functions Gk and Tk 
coincide, see for instance [44] and the references there. However, a rigorous formu- 
lation as a minimal problem for the free energy reveals that Qk — F'k * s appropriate. 
This topic has been thoroughly investigated for analogous phase separation prob- 
lems, see [40, 41, 22, 23], see also [18] and [24]. In order to cover both cases we 
regard independent functions Gk an d Tk- 

2.3 Assumption. Mathematically, we demand that the distribution functions 
Gk, k = 1,2, are defined on the real line, take positive values, and are either exponen- 
tials, or twice continuously differentiable and polynomially bounded. Moreover, J-~[, 
T' 2 are strictly positive on R. In the sequel we will call such distribution functions 
'admissible.' This includes Boltzmann statistics, as well as Fermi-Dirac statistics 
(see (2.8)). 

Let us comment on the (effective) band edges bk and the (effective) densities of 
states pk, see (2.5) and (2.6): Basically the band edge offsets bk and the effective 
band edge densities of states pk are material parameters. In a heterogeneous semi- 
conductor device they are generically piecewise constant on the spatial domain Q. 
As Assumption 3.7 reveals, we cannot cope with such a situation as far as the band 
edges bk are concerned. However, in the case of Boltzmann statistics one can rewrite 
(2.5) and (2.6) as 



with modified effective densities of states and identically vanishing band edge offsets. 
In the case of Fermi-Dirac statistics this reformulation is not possible and one has to 
recourse to some approximation of the bk by functions confirming to Assumption 3.7. 
Discontinuities of the band edge offsets up to now seem to be an obstacle in whatever 
approach to solutions of van Roosbroeck's equations, if the statistical distribution 
function is not an exponential, see for instance [19]. 

There are compound multiscale models of semiconductor devices such that the ef- 
fective band edges and the effective densities of states result by upscaling from 
quantum mechanical models for the electronic structure in heterogeneous semicon- 
ductor materials, see [2, 3, 35]. In view of an offline coupling to electronic structure 
calculations we allow for an explicit dependence of pk, and bk on time. 

2.4 Reaction rates 

The reaction terms on the right hand side of the current-continuity equations can 
be rather general functions of time, of the electrostatic potential, and of the vector 
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of the electrochemical potentials. r\ and r 2 describes the production of holes and 
electrons, respectively — generation or annihilation, depending on the sign of the 
reaction term. Usually van Roosbroeck's system comprises only recombination of 
electrons and holes: r = r\ = r 2 . We have formulated the equations in a more general 
way, in order to include also coupling terms to other equations of a superordinated 
compound model. That is why we also allow for an explicit time dependency of the 
reaction rates. 

Our formulation of the reaction rates, in particular, includes a variety of mod- 
els for the recombination and generation of electrons-hole pairs in semiconductors. 
This covers non-radiative recombination of electrons and holes like the Shockley- 
Read-Hall recombination due to phonon transition and Auger recombination. But, 
radiative recombination (photon transition), both spontaneous and stimulated, is 
also included. Mathematical models for stimulated optical recombination typically 
require the solution of additional equations for the optical field. Thus, the recombi- 
nation rate may be a non-local operator. Moreover, by coupling van-Roosbroecks 
system to the optical field some additional control of this optical field may also 
interact with the internal electronics. For instance, in modeling and simulation of 
edge-emitting multiple-quantum-well lasers van-Roosbroeck's system augmented 
by some Helmholtz equation often serves as a transversal (to the light beam) model, 
and a control of the optical field is exercised by a master equation or some model 
for the longitudinal (on the axis of the light beam) behavior of the laser, see for 
instance [51, 2, 3]. 

Modeling recombination of electron-hole pairs in semiconductor material is an art 
in itself, see for instance [36]. However, for illustration, let us list some common 
recombination models, see for instance [44, 11] and the references cited there. 

Shockley-Read-Hall recombination (phonon transitions): 



where rij is the intrinsic carrier density, n 1; n 2 are reference densities, and T\, t 2 are 
the lifetimes of holes and electrons, respectively, ni, n±, n 2 , and n, r 2 are parameters 
of the semiconductor material; thus, depend on the space variable, and ultimately, 
also on time. 

Auger recombination (three particle transitions): 



where c 1 uger and c 2 uger are the Auger capture coefficients of holes and electrons, 
respectively, in the semiconductor material. 
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Stimulated optical recombination: 

m 2 



n = r 2 = r stim = f( a S 



"Ml 2 ' 

where / additionally depends on the vector of the densities, and on the vector of the 
electrochemical potentials, aj, ipj are the eigenpairs of a scalar Helmholtz-operator: 

Atpj + e(ui, u 2 )i>j = o-jipj. 

In laser modeling each eigenpair corresponds to an optical (TE) mode of the laser 
and \ipj\ 2 is the intensity of the electrical field of the Oj— mode, e is the dielectric 
permittivity (for the optical field); it depends on the density of electrons and holes. 
The scalar Helmholtz-equation originates from the Maxwell equations for the optical 
field [50]. 

The functional analytic requirements on the reaction terms will be established in 
Assumption 3.6. 



3 Mathematical prerequisites 

In this section we introduce some mathematical terminology and make precise as- 
sumptions about the problem. 



3.1 General Assumptions 

For a Banach space X we denote its norm by and the value of a bounded linear 
functional ip* on X in if) G X by (ip* \ip)x- If X is a Hilbert space, identified with 
its dual, then (• | -)x is the scalar product in X. Just in case X is the space M 2 , the 
scalar product of a, b G R 2 is written as a -b. Upright X denotes the direct sum X®X 
of slanted X with itself. B(X; Y) is the space of linear, bounded operators from X 
into Y, where X and Y are Banach spaces. We abbreviate B(X) = B(X; X) and 
we denote by Boo(X) the space of linear, compact operators on the Banach space 
X. The notation [X, Y]g means the complex interpolation space of X and Y to the 
index 9 G [0,1]. The (distributional) V-calculus applies. If ^ is a (differentiable) 
function on an interval taking its values in a Banach space, then if)' always indicates 
its derivative. 



3.2 Spatial Domains 

Throughout this paper we assume that f2 as well as f2 are bounded Lipschitz domains 
in M 2 , see [25, Ch. 1]. By ' we denote the operator which extends any function 
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defined on Q by zero to a function defined on f2. Conversely, j_ denotes the operator 
which restricts any function defined on Q to Q. The operators T and j are adjoint to 
each other with respect to the duality induced by the usual scalar product in spaces 
of square integrable functions. 

With respect to the marked out Neumann boundary parts T C d£l and T C <9f2 of 
the boundary of Q and Q we assume each being the union of a finite set of open 
arc pieces such that no connected component of dQ \ T and dQ \ T consists only 
of a single point. We denote the parts of the boundary where Dirichlet boundary 
conditions are imposed by D = interior(<9f2 \ T) and D = interior (dQ \ T). 

3.3 Function spaces and linear elliptic operators 

We exemplarily define spaces of real- valued functions on spatial domains with respect 
to the bounded domain Q C M 2 and its boundary. Spaces of functions on Q and 
parts of its boundary may be similarly defined and are denoted by hatted symbols. 

If r G [1, oof, then If is the space of real, Lebesgue measurable, r-integrable functions 
on Q and L°° is the space of real, Lebesgue measurable, essentially bounded functions 
on f2. W 1,r is the usual Sobolev space W 1,r ({l), see for instance [46]. W^ r is the 
closure in W 1,r of 

{^|n : ^6QK 2 ), supp ^ n (dQ \ T) = 0} , 

i.e. Wy' t consists of all functions from W 1,r with vanishing trace on D. ' r 
denotes the dual of W^' r , where 1/r + l/r' = 1. (■ | -) w h2 is the dual pairing 

1 2 1 2 ^ 

between W T ' and W T ' . Correspondingly, the divergence for a vector of square 
integrable functions is defined in the following way: If j G L 2 , then V • j G ' is 
given by 

(V-j\1>) w ^ = - /j-V^dx, ^eW^ 2 . (3.1) 
r Jn 

a is the natural arc measure on the boundary of Q. We denote by L°°(dQ) and 
L r (dQ), the spaces of cr-measurable, essentially bounded, and r-integrable, r G 
[1, oof, functions on dQ, respectively. Moreover, W s ' r (dQ) denotes the Sobolev space 
of fractional order s G]0, 1] and integrability exponent r G [l,oo[ on dfl, see [25, 
Ch. 1]. Mutatis mutandis for functions on a-measurable, relatively open parts of 

on. 

Let us now define in a strict sense the (linear) Poisson operator and the elliptic 
operators governing the current continuity equations. 

3.1 Definition. We define the Poisson operator —V ■ eV : W 1,2 — » W~ 1,2 by 

(-V ■ eV^i | ^ 2 )w^ = I ■ V^ 2 dx+ I e f ^ 2 da, (3.2) 

r Jn Jf 



Preprint 1189, Weierstrass Institute for Applied Analysis and Stochastics, Berlin 2006 



Classical solutions of drift-diffusion equations 



11 



' — "i 9 " — ^1 2 " — ^1 2 

for i/ji £ W ' and tp2 £ W~' . Vq denotes the restriction of —V • eV to W~' ; we 
denote the maximal restriction of Vo to any range space which continuously embeds 
into ' by the same symbol TV 

3.2 Definition. With respect to a function q £ L°° we define the operators 

- V ■ <^V : W 1 ' 2 -> Wf 1 ' 3 , fc = 1, 2, by 

(-V ■ ^/i fc V^! I = / ? /i fe V^i • V^ 2 dz, e W 1,a , ^2 e W?' 2 . 

r in 

If, in particular, ? = 1, then we simply write d k for —V ■ /UfcV. Moreover, we denote 

i 2 12 12 

the restriction of a& to the space W T ' by a^, i.e. a k : W r ' — > W r ' . 

3.3 Proposition, (see [26] and [27]) There is a number q > 2 (depending on Q, e 
and Y) such that for all q £ [2,g] the operator Vq : W~ q — > W~ 1,q is a topological 
isomorphism. Moreover, there is a q > 2 (depending on Q, \xi and T) such that 
for all q £ [2, q] the operators a k : W^ q — > ' q provide topological isomorphisms, 
and additionally, generate analytic semigroups on ,q - 

3.4 Definition. From now on we fix a number q £]2, min(4, q, q) [ and define p = |. 
With respect to this p we define the operators 

A k : V a fc ^, ^ e £>* = dom(A fe ) = {V> £ W^ 2 : a fc ^ £ L p } , = 1, 2, 
A : D -> A = £) , 2? = dom(A) = £>x © P 2 L" . 

3.5 Remark. If ip £ X>j., fc = 1, 2, then v-(//feV , 0)|r = in the sense of distributions, 
see for instance [5, Ch. 1.2] or [16, Ch.1.2]. 

After having fixed the number q and, correspondingly, the space LP ', we will now 
formulate our mathematical requirements on the reaction terms: 

3.6 Assumption. The reaction terms r k , fc = 1,2, are mappings 

r k : [T ,Ti] x W 1,q x W 1 ' 9 -> L p . 

Moreover, we assume that there is a real number rj £]0, 1] and for any bounded 
subset M C W l ' q © W 1 ' 9 a constant r A/ such that 

\\r k (t,v,ip) -r k (i,v,ip)\\ LP 

< r M (\t - i\ v + \\v - v\\^ ltq + - ^Hw 1 ^) , 

t,ie [To,Tx], (^), (vJ)eM. 

3.7 Assumption. The functions b k : [T ,Ti] — > W 1,q , k = 1,2, are Holder contin- 
uous. Moreover, they are Holder continuously differentiable when considered as L p 
valued. 



Preprint 1189, Weierstrass Institute for Applied Analysis and Stochastics, Berlin 2006 



12 



H.-Chr. Kaiser, H. Neidhardt, J. Rehberg 



3.4 Representation of Dirichlet boundary values 

For setting up the Poisson and current-continuity equations in appropriate function 
spaces we must split up the solution into parts, where one part represents the inho- 
mogeneous Dirichlet boundary values and (f>r>,k, k — 1,2. In this section we treat 
of just this representation. We make the following assumptions about the Dirichlet 
boundary values of the electrochemical potentials <pk, k — 1,2, and for their initial 
values, see (2.3), (2.4). 

3.8 Assumption. There is a Holder continuous function 

$ = ($ 1 ,$ 2 ):[T ,T 1 ]^W 1 '", A; = 1,2, 
such that for all t e [To, Xi] 

a k $ k (t) = (3.3) 
tr($ fc (i))| D =<l> D , k (t) (3.4) 

Moreover, we assume, that each k — 1, 2, — as a function with values in LP - 
is differentiable and its derivative is Holder continuous. 

3.9 Remark. It should be noted that (3.3) and the definition of the operators 
imply v ■ /ifcV$fc = on T in the distributional sense, see for instance [5, Ch. 1.2] or 
[16, Ch. II. 2]. This implies for the current densities (2.9) that v ■ = on T in the 
distributional sense, provided that \k £ W 1,q . 

We will now give a sufficient condition on tp^^ for the existence of a with the 
assumed properties. 

3.10 Lemma. 1. If ip e W 1-1 / M (D), then there is a unique function e W ' q 
fulfilling 

a k y = 0, and tr(^) | = -0. 

2. Ifip : [Tq, Tj] — > H^ 1 ~ 1 /' 3 ' <7 (D) is Holder continuous with index r], then the function 
\1/ : [Tq,Ti] — > W 1,q which is given for each t G [T ,Ti] by item 1 is also Holder 
continuous with index rj. Moreover, if tp — as a function with values in W l l 2,2 (D) 
- is Holder continuously differentiable with Holder index r\, then is Holder con- 
tinuously differentiable with Holder index rj. 

Proof. Let ex : W 1 ~ l ^ q,q {Y)) — > W 1 ~ l ^ q,q {dVt) be a linear and continuous extension 
operator, and let tr _1 be a linear and continuous right inverse of the trace operator 
tr : W 1,q {&) -> W 1 ~ 1 / q ' q (dn). Such operators exist according to [25, Thm 1.4.3.1] 
and [25, Thm 1.5.1.3], respectively. Thus, tr -1 oexip e W 1,q . Moreover, let ip be the 
solution of the differential equation 

akip = dk o tr -1 oextjj (3.5) 
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in W^ 9 . This solution exists and is unique because the right hand side of (3.5) is 
from Wy ,q and the operators a k are isomorphisms from W^ 9 onto ' q . We now 
define 

* = tr^oex^-^- (3.6) 
The asserted properties of \I/ follow directly from the construction. 

The second assertion is proved by observing that all steps in the first part of the 
proof depend linearly on the datum. □ 

3.11 Assumption. We assume that the initial values $° belong to W 1,q , k = 1,2. 
Moreover, there is a # G]l/2 + 1/g, 1[ such that for each of the initial values the 
difference $° — $fc(T ) belongs to the complex interpolation space [L p ,V k ] e . 

3.12 Remark. For all 9 G]l/2 + 1/q, 1[ the space [L p ,Dk]g compactly embeds into 
W}> q ^ L°°, see [31, Thm. 5.2]. 

With respect to the inhomogeneous terms y2g and (pf- in the boundary conditions of 
Poisson's equation (2.1) we make the following assumptions. 

3.13 Assumption. There is a Holder continuous function (p Q : [T ,Ti] — > W 1,q such 
that f — as a function from [To, Ti] into L p — is Holder continuously differentiable. 
For all t G [T , Ti] it holds true 

-V-evVo(0 = 0, (3.7) 
tr(^ (t))| fi = ^g(t). (3.8) 

The function 

[T ,T 1 ]9t^^ F (t)eL oo (f) 
is differentiable and possesses a Holder continuous derivative. 

3.14 Remark. Similar to Lemma 3.10 it is possible to give a sufficient condition on 
the existence of a representing function t i— > (p Q {t) which only rests on the function 
t i — > <PQ{t). We do not carry out this here. 

3.15 Remark. For all t G [To,Ti] we extend Pf{t) by zero to a a-measurable, 
essentially bounded function on dQ. Due to the continuous embedding 

Wh q ' ^ W u ^ W 1 ~ 1 t q ' ,q ' (dti) ^ L q '(dh), 

see [25, Thm 1.5.1.3], there is a continuous embedding 

L°°(dh) ^ L q (dQ) ^ W~ l ' q . 

Thus, y?f(t), * e Po>Ti] can be regarded as an element of W~ 1,q . We denote iff- as 
a function from [T ,Ti] into W- 1 ' 9 by ip 9 . The Holder continuous differentiability of 

<Pf- entails the Holder continuous differentiability of <p, : [To,Ti] — > W~ ' q with the 
same Holder exponent. 
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3.5 The linear Poisson equation 

Let us assume the following about d — the doping profile (or control parameter) on 
the right hand side of Poisson's equation (2.1). 

3.16 Assumption. The function d : [T ,Xi] — > W~ l,q is continuously different iable 
with Holder continuous derivative. We define a "generalized doping" 

d : [T , T x ] - WT 1 ' 9 by d(t) = d(t) +<p.{t), t G [T , T,]. (3.9) 
We now define what is a solution of Poisson's equation (2.1). 

3.17 Definition. Let Uf- G ' q , k — 1, 2 be given. We say that ip is a solution of 
Poisson's equation (2.1) at t G [T ,Ti], if 

fi=ip + ip {t), (3.10) 

and tp G Hp' 9 is the unique solution of 

V <p = d(t) + Ul - u 2 . (3.11) 

ip and <p depend parametrically on t, U\, and Ui- If convenient, we indicate the 
dependence on t by writing ip(t) and <p(t), respectively. 

3.18 Remark. With respect to the boundary conditions in (2.1) it should be noted 
that (3.8) and the property ip G W~ q give = ipQ. Additionally, if d, u\, and 
u 2 belong to the space L 1 , then (3.9), (3.10) and (3.11) together with (3.7) imply 
v ■ [eSJp] + Ef <p = <pf(t), see for instance [5, Ch. 1.2] or [16, Ch. II. 2]. 

Throughout this section we demand several times Holder continuity of functions 
and/or their derivatives. Clearly, there is a common Holder exponent which we will 
denote from now on by 77. 

4 Precise Formulation of the Problem 

We are now going to define the problem outlined in §2. 

4.1 Definition. We say the van Roosbroeck system admits a local in time solution, 
if there is a time T G]T , T x ] and (!p, 0) = (!p, l5 2 ) such that 

f (T ) = (0i(T o ), MT )) = ($?, $°) G W 1 ' 9 , (4.1) 

p ^ p - Vo e C([T , T]-Wh q ) n C\]T , T[; H/i' 9 ) (4.2) 
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= — $ G C\]T , T[, I/) n CQTo, T], P) n C([T , T], [I/, 2?],), (4.3) 
fulfill the Poisson equation and the current continuity equations: 

V {<p(t)) = d(t) + V(t) - T w 2 (t) t e [T , T], (4.4) 

<(*)- V-i fc (O=r fc (t,^),0(t)), fc = l,2, te]T ,T[. (4.5) 
The carrier densities and the current densities are given by 

u k (t)=p k (t)T k ( X k(t)), (4.6) 

jk(t) = Qk(Xk(t))(M k VMt)> ( 4 - 7 ) 
X*W = ?*(*) + (-l)*iTO + &*(*)■ 



and satisfy 



u k e C([T ,T],L°°) n C^lTo,^,^), (4-9) 

iteC^nn (4.10) 

V-j fe eC(]To,T],I7) (4.11) 



for fc = 1,2. 



5 Reformulation as a quasi-linear parabolic sys- 
tem 

In this section we provide the tools to rewrite the problem from Definition 4.1 as 
a quasi-linear system for the continuity equations. To that end we eliminate the 
electrostatic potential from the continuity equations. Replacing the carrier densities 
u\ and ui on the right hand side of (4.4) by (4.6) making use of (4.8) and (3.10) one 
obtains a nonlinear Poisson equation for (p. We solve this equation with respect to 
prescribed parameters b k and (f) k , k = 1,2, which we will assume here to be from L°°. 
This way to decouple van Roosbroeck's equations into a nonlinear Poisson equation 
and a system of parabolic equations is also one of the fundamental approaches to 
the numerical solution of the van Roosbroeck system. It is due to Gummel [28] and 
was the first reliable numerical technique to solve these equations for carriers in an 
operating semiconductor device structure. 



5.1 The nonlinear Poisson equation 

We are now going to prove the unique solvability of the nonlinear Poisson equation 
and some properties of its solution. First we show that the supposed admissibility 
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of the carrier distribution functions T\. ensures that the relation between a potential 
and its corresponding carrier density is monotone and even continuously differen- 
tiable when considered between adequate spaces. 

5.1 Lemma. Let p and g be from L°° and T = Tk be an admissible carrier distri- 
bution function, see Assumption 2.3. 

1. The operator 

Wh 2 3h^ P F{ 9 + i h)eL 2 (5.1) 



is well defined, continuous and bounded. Its composition with the embedding L 2 
W~ ' is monotone. 

2. The Nemyckii operator 



L°°3h^pF(g + l h) 

induced by the function 

fi x R 3 (x, s) i — > p(x)F(g(x) + s), 

maps L°° continuously into itself and is even continuously differentiable. Its Frechet 
derivative at h e L°° is the multiplication operator given by the essentially bounded 
function 

fl9ii — > p{x) F'(g(x) + h(x)). (5.2) 

Proof. Indeed, the assumption that the carrier distribution functions should be ad- 
missible assures that the operator (5.1) is well defined, continuous and bounded, see 
[47] for the case of an exponential, and see [1, Chapter 3] for the case of a polyno- 
mially bounded function. The asserted monotonicity follows from the monotonicity 
of the function T and the fact that the duality between W~' and W~ ' is the 
extension of the L? duality: 

r 

= / ^ P F{g + - ^pT{g + x h 2 )) (h - h 2 ) dx 
Jn 

= I {pF{9 + ihi) - pT{g + { h 2 )) - t h 2 ) dx > for all h u h 2 G W± 2 . 
Jn 

The second assertion follows from a result by Groger and Recke, see [42, Thm 5.1]. 

□ 

5.2 Corollary. The mapping 

W 1 «3h^lpf(g + l h) 
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takes its values in L°° and is also continuously differentiable. Its derivative at a 
point h G W 1 ' 9 equals the multiplication operator which is induced by the function 
ipFig + M. 

5.3 Theorem. Under Assumption 2.3 on the distribution functions T\, T 2 and 
Assumption 2.1 the following statements are true: 

1. For any pair of functions z — (z\, z 2 ) G L°° the operator 

if 1 — ► V (f - Vi-^i^i - i<p) + ^P2^2{z 2 + i<p) (5.3) 

is strongly monotone and continuous from W~ 2 to W~ 1,2 , where the operator Vq is 
according to Definition 3.1. The monotonicity constant of (5.3) is a least that of 
V . 

2. For all f G W^ 1,2 and z = (zi, z 2 ) G L°° the nonlinear Poisson equation 

V if - T pi-^i(^i - if) + ^P2^2{z 2 + i<p) = f (5.4) 

admits exactly one solution <p which we denote by C(f,z). This solution belongs to 
W~' and satisfies the estimate 

\\<p\\ww < — HW^i) - ^P%T 2 {z 2 ) + fWwz 1 - 2 ' 
where m is the monotonicity constant ofVo- 

3. The maximal restriction of the operator (5.3) to the range space Wf 1,q has the 
domain W~' q . Moreover, if M is a bounded subset of W~ 1,q © L°°, then the set 
{£(f,z) : (f,z) G M} is bounded m W~' q . 

4- The mapping L : W~ ,q ' ©L°° — > W~ q is continuously differentiable. Let (F, Z) = 
(F, Zi, Z 2 ) be from W~ 1,q © L°°; we define the function 

M k d = f WftS* + (-1)\C(F, Z)), (5.5) 

and we also denote the corresponding multiplication operator on Q by Af k . Then the 
Frechet derivative dC at a point (F, Z) = (F, Zi, Z 2 ) is the bounded linear mapping 
given by 

[dC{F, Z)\ (/, z) = {V + M + ^y 1 (f + M T *i - M 2 h 2 ) , k = 1, 2 (5.6) 
for all (f, z) = (f, ( Zl , z 2 )) G W f hq © L°° . 

5. The norm of dC(F, Z) G B{Wp q © L°°; Wp q ) can be estimated as follows: 

\\dC(F, Z)|| B ^-l,g eL oo.^l,9) 



< 2 ll^ 1 IU(L 2; ^)V / IIM+Ar 2 || i »||AA 1 +AA 2 || il + IIPq" 1 ! 



B(W~ 1 ' q ;W~' q ) 



+ ll^o 1 \\b(Z 2 -w± 9 )V\\Ni + N 2 \\l°°\\Vq 1 2 || s( ^ 



r 1,9 ii a ) 
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Proof. 1. The assumption that D is not empty or e-p is positive on a set of positive 

arc measure ensures that the operator P is strongly monotone. Thus, taking into 

account Lemma 5.1, the mapping (5.3) is strongly monotone and continuous from 

W~' 2 to W- 1 ' 2 . 
r r 

2. The second assertion follows from the first one by standard results on monotone 
operators, see for instance [16]. 

3. For / G W~ ' the solution £(/, z) is from W~' and hence, 

- Vi^i (zi - iC(f, z)) + V 2 jr 2 (* 2 + iiC (/, z )) eL 2 ^ W^ q , 

see Lemma 5.1. By the second assertion of the theorem, the set 

{£(/, z) : (/, ^)eM} is bounded in W~> 2 . 

From this we conclude again by Lemma 5.1 that the set 

{^p 1 ^ 1 (z 1 - l /:{f,z))-^ 2 ^(z 2 + l /:(f,z)) : (f,z)EM} 

is bounded in L 2 , and hence, is bounded in W~ 1,q . Thus, the set 

{^p 1 ^ 1 ( Zl - l C(f,z))-^p 2 J' 2 (z 2 + l C(f,z))+f : (f,z)EM) 

is also bounded in W~ 1,q . Consequently, the image of this set under Vq 1 is bounded 
in W~' q . 

4. We define an auxiliary mapping K : W~ ,q © W~ 1,q © L°° -> 1,9 by 

£(<P, /, Z ) = ^OV - T pi-^l(^l - 1<P) + ] P^{z 2 + l<p) ~ f 

such that )C(jC(f,z),f,z) = for all / G WT 1 ' 9 and all z G L°°. The assertion 
follows from the Implicit Function Theorem if we can prove that /C is continuously 
differentiable and the partial derivative with respect to tp is a topological isomor- 
phism between W~' q and WT l ' q . For any tp G Wi' 9 , / G WT 1 ' 9 , and z G L°° the 
partial derivatives of JC are given by 

2 

<yC(^,/,z) = V + ^PkH(zk + (-l)\v)eB(W± q ;W~ 1 > q ), (5.7) 

fe=i 

d f K{ipJ,z) = -IeB(W^' q -,W^' q ), (5.8) 
dJCfaM = (-l) k Tp k Fl(z k + (-l)\tp)EL™^B(L™-W~ 1 ' q ) (5.9) 

and they are continuous, see Lemma 5.1 and [42, §5]. 
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Now we consider the equation 

2 



(5.10) 



k=l 



Because Ylk=i Vfc-^fc(^fc + ( — ~^) k i^) is a positive function from L°°, (5.10) has exactly 
one solution ip e W~' by the Lax- Milgram- Lemma. Moreover, 



1,9 



fc=l 



and P : W- 9 -> WT 1 ' 9 is a topological isomorphism. Thus, a rearrangement of 
terms in (5.10) gives ip G W~' q . 

5. We now estimate the Frechet derivative (5.6): 

r 

< \\{V Q +N 1 +N 2 )- 1 f\\^ 

r 

+ \\{V Q +N 1 +N 2 )- l {Ujz l -Njz 2 )\\^ q . (5.11) 

r 

We treat the right hand side terms separately; for the second addend one obtains 



r 



< 



. 9 ^ l9 Nlra, (5-12) 



S(L 2 ;W~' 9 ) 



where the function g £ L 2 is defined by 

dcf M(a;)zi(x) - M 2 {x)z 2 (x) 



9{x) 



y/Afi(x)+Af 2 (x) 



for i6Q. 



(5.13) 



Please note that the functions A4 are strictly positive almost everywhere in Q due 
to the positivity of the distribution functions and Assumption 2.1. For the function 
g in (5.13) one has the following bound: 



9\\v < \ ||M +M 2 \\ti (IkilU- + 



Making use of the operator identity 



(Vo + Mi + M 2 ) ~ l = TV - Pc i (M + (Po 



-i 



(5.14) 
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one obtains 



B(L 2 ;W~' q ) 



< 



S(L 2 ;VKi' 9 ) 



+ 



< ll^ 1 



B(L 2 ;Wi' ? ) 



lB(L 2 ;W^) 



X 1 + 



2 

B(£ 2 ) 



We note that 



1/2 



B(Z 2 ) 



< 1 



(5.15) 



because the bounded multiplication operator A/i + A/i? is form subordinated to Po + 
A/i + A/2, see for instance [34, VI. 2. 6]. Thus, we get for the second addend of (5.11): 

HCPo+M + M) - W*i - A/^) lb: 



r 



- 2 \\V ( 



'X^M") yJ\\^+^\\z«>yJWi+^2hi (INU~ + INU~) (5.16) 



Applying (5.14) to the first term on the right hand side of (5.11) we find 



II'- 



-1,9 



+ ll^o ^Ib^jwM) H(M + M)(^o + M + A/" 2 ) 'H^-M.^ 



r 



(5.17) 



The terms |P n 1 1 



and 1 1 1 e(Z 2 ■ vk 1 ' 9 ) are ^ n ^ e - As f° r the remaining 



term 

IKm+a/sXPo+a/^+a^)- 1 ! 



< ^/||A/i+A4|| £o c VKTM(v + m + a/" 2 )- 1/2 

(Po+A^+A^)" 1 ^ 



B(L 2 ) 



B(L 2 ) 



-1/2 



B(WT 1,9 ;L 2 ) 



we note that 



.1/2 



1.2 



is finite, since W~ ,q embeds continuously into W~ ' 

-v^- 1 "^) ' r * r 

and P 1/2 : L 2 -> W-T 1,2 is a topological isomorphism. Again, P is form subordinated 
to "Po + A/i + A/2. Hence, besides (5.15) one has 

B(L 2 ) - 1- 

Thus, we get from (5.17): 
\\(Vo + M l +M 2 r 1 f\\ m , q <\\Vo 1 



\B{WI 1 ' q ;W~ q ) 



r 



' 



B(2 2 ;W / i' 9 ) 



IIM+jsr 2 |h 



-1/2 



B(Wp 1,? ;£ 2 ) 



r 



(5.1* 
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Inserting (5.16) and (5.18) into (5.11) finishes the proof. □ 

5.4 Corollary. Let the assumptions of Theorem 5.3 be satisfied. Then holds true: 

1. The mapping C : W~ ,q © L°° — > W~ q is boundedly Lipschitzian, i.e. for any 
bounded subset M C W~ x ' q © L°° there is a constant Cm such that 

||£C/» - £Uiz)\\ w i, q < £m ( ||/ - fWwz 1 " + II 2 ~~ 



for all{f,z), (f,z)eM. 

2. Let additionally Assumption 3.16 be satisfied. If 

z={z u z 2 ) G C([T , T], L°°) n C\]T Q , T[, L p ), 



toera toe function [T ,T] 9 t i-> y>(t) G Wp' 9 (/wen 6y ^(t) = C(d(t),z(t)) G 
continuous, and continuously differentiable on ]Tq,T[. Its derivative is 

<f/(t) = [d£(d{t),z{t))] (d\t),z'(t)) 

= (Vo + Afi + A^)" 1 (d'(t) + MM - A/" 2 T 4) , 

where Mk is again defined by (5.5) — there (F, Z) specified as (d(t), z(t)) . 
5.2 Derivation of the quasi- linear system 

We start now with the reformulation of the van Roosbroeck system as defined in 
Definition 4.1 as a quasi-linear parabolic system for the continuity equations. The 
aim of eliminating the electrostatic potential in mind, we first look for a substitute 
for its time derivative. In order to achieve this, we formally differentiate Poisson's 
equation (4.4) with respect to time. This gives 

TV = d' + \u' l - u 2 ). (5.19) 

From (4.5) one obtains 

u[ - u' 2 = V • ji - V • j 2 + n(i, <p, 4>) - r 2 (t, <p, 4>). (5.20) 

Inserting (5.20) into (5.19), one gets 

TV = df + T (V • jx - V • j 2 + n(t, <p, 0) - r 2 (t, <p, 4>)) . (5.21) 

Just in case, r = r\ = r 2 is only recombination, this is precisely the well known 
conservation law for the total current, see [11]. Clearly, (5.21) leads to 

If' = iK 1 (d' + T (V • 3l - V ■ j 2 + n(t, <p, 4>) - r 2 (t, tp, 0))) . (5.22) 



is 
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Now we differentiate (4.6) (with (4.8)) with respect to time and obtain 

< = Pkf&k + (-l) fc i£ + 6*) R + (-1) V + b' k ] 

+ p' k FSk + (-l) k & + h), k = l,2, (5.23) 

Pending further notice we do not write out the argument (p k + (— l) k iP + b k of the 
distribution function T k and its derivative. We also abstain from drawing out the 
argument of the reaction terms r k . According to (3.10) we split ip' = tp' + <p' and 
insert (5.23) into the current continuity equation (4.5). Thus, we find 

+ {-l)\ip']p k r' k -V- Jk = r k - [(-1) Vo + &'J p k T' k - p' k F k , k = 1, 2. 
Using (5.22) we get further 

- V • j fc + (-1) W^o" 1 K + T (V • ji - V • h + n - r 2 )) 

= r„- [{-l)\<p' + b' k ]p k F' h -p' k F k , fc = l,2. 

Dividing this by Pfc^. we obtain 

/ i + {p^p iPx -{P^r 2P2 
\ -iP^r lPl i + {P^r 2 p 2 





This evolution equation can be written in the condensed form 

0' - [/ + Z(t, $)]E{t, 0) V • j = Y(t, 0) (5.24) 

where = (0i,02) and V • j = (V • ji, V ■ j 2 ). Moreover, / denotes the identity. 
The coefficients Z, E, and Y are given in the following way: First we split off the 
Dirichlet inhomogeneities of tp in the sense of §3.4 and we replace <p by the solution 
of the nonlinear Poisson equation, see Theorem 5.3. With respect to an arbitrary 
ip = (ip!, ip 2 ) <E W 1 ' 9 we set 

Qk{t,ip) = 4> k + (-l)V(rf(t),^(t)) + (-l)W*) +h(t), k= 1,2, (5.25) 

where z = (zi, z 2 ) with 

z k (t) =i> k + {-l)\p {t) + b k {t), = 1,2. (5.26) 

Now we define 

7(t ,M d = f ( I^O _1T ^(Ql(i^))PlW -|V T ^2(Q2M))P2(*) \ 27 x 



. 1( n^Ktp^ + ^jt),^) 
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and finally 

r(*,V)= [l + Z{t^)]E{t^)R{t^)-X{t^), (5.30) 
where X(t,ip) = (X 1 (t,ip),X 2 (t,rl>)) with 

x k {t^) d = 4 (-l)*^ 1 ^*) + */„(*)) + b' k (t) + (5 31) 

k = 1,2. Please note 

ZM)E{t^)=( (5.32) 

Next we apply the definition (2.9) of the currents j k and get 

V • j k = V ■ (£ fc (0 fc + (-1) V + (-l) fe i^o + b k )pL k V$ k ), k = 1, 2, 
or in shorter notation 

V-j = V-G(t,0)/A70, (5.33) 
where — see also (5.25) and (2.9) — 

G(t,1>) d = f ( Gl( ^ G2{ t ^) , G k {t^) d ^g k (Q k (t^)). (5.34) 

Now, putting together (5.33) and (5.24) we obtain in conclusion the evolution equa- 
tion 

?'-[/ + £(t, 4>)] E(t, 4>) V ■ G(t, 0)a*V0 = F(t, 0) (5.35) 

which has to be complemented by the boundary conditions (2.3) and the initial 
condition (2.4), see also Remark 3.9. 



6 The quasi-linear parabolic equation 



Evolution equations of the type (5.35) were investigated in [31]: (5.35) has a unique, 
local in time solution, if the functions E, G, Z and Y defined by (5.28), (5.34), (5.27) 
and (5.30), respectively, satisfy the following conditions. 

6.1 Assumption. With respect to q e]2,oo[ and p = q/2, as specified in Defi- 
nition 3.4, there is an 77 g]0, 1] and further for any bounded set M C W 1 ' 9 exist 
positive constants E M , Gm, Y m , and Z M such that the mappings 



E 
G 
Z 
Y 



[T 0j Fx] x W 1 ' 9 
[To, Ti] x W 1 ' 9 



L°°, 
W 1 ' 9 , 

17 



(6.1) 
(6.2) 
(6.3) 
(6.4) 
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satisfy the conditions 

min inf vraimin ib)(x) > (6.5) 

fc=l,2 ie[T ,Ti] rrGf! V 7 V ^ 

min inf vraimin il})(x) > (6.6) 
fc=i,2 te[T ,Ti] xen 

and for all t, t G [T , Ti] and all ip,ip E M: 

\\E(t,i/;)-E(i,<4>)\\ L °o < E M (\t-i\^ + \\^-t[j\\ w ^, q ), (6.7) 

||G(t,^)-G(t,^)|| wl , 9 < G M (|t-^ + ||^-^|| wl „), (6.8) 

||Z(t,^)-Z(t,^)|| B(LP) < Z M (|t-t|"+||^-^|| wl , ? ), (6.9) 

\\Y(t,i;)-Y(i,4>)\\ LP < Y M (\t-i\*> + \\i,-4>\\ w i, q ). (6.10) 

6.2 Definition. Let the Assumptions 3.8 and 6.1 be satisfied. Further, let A : 
T> — > L p be the operator from Definition 3.4 and let V be a Banach space such 
that T> ^ V ^ W 1 ' 9 . We say the evolution equation (5.35) with initial condition 
4>{Tq) = $° G W 1 ' 9 has a unique local solution = + $ with respect to V if 
$° — $(T ) G V implies the existence of a number T g]Tq,Ti] such that the initial 
value problem 

<j>'{t) + [I + Z (t, <j>{t) +<S>(t))]E(t,<f> + $(t))G (t, 0(t) +$(*)) A 0(t) 

= F(t,0(t) + $(£)) -&(t) + J(t,<f>(t)), 0(T o ) = $°-$(T o ) (6.11) 

admits a unique solution 

G C^QTo, T[, I/) n C(]T , T], P) n C([T , T],V). (6.12) 

For (t,-0) G [T ,Ti] x Wp' 9 the term J in (6.11) is given by 

j(t,if>) = [i + + + $(*)) vc(t,v + $(*)) -/iV(^ + $(t)). 

6.3 Remark. We have to clarify the relation between (5.35) and (6.11). If = 0+$ 
is a solution in the sense of Definition 6.2, then 

V-G(t,0)/iV0 = G(t,0)A0 + VG(t,0) -/iV0 (6.13) 

is satisfied, which allows to rewrite (6.11) in the form (5.35). 

6.4 Remark. If = (0i,02) is a solution of (5.35) in the sense of Definition 6.2, 
then 

tr(0 fc (t))| D = tr($ fc (t))| D = Difc (t), A; = 1,2, t G [T , T]. 
The Neumann boundary condition 

= v-VikVMt)\ T = u-ti k V$ k (t)\ T , k = 1,2, te[T ,T], 

holds in the distributional sense, see Remark 3.9. 
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6.5 Proposition. (See [31].) Let the Assumptions 3.8 and 6.1 be satisfied. For 
each 7 £ ]~ + ^,l[ the initial value problem (5.35) with initial value $° £ W ' q 

has a unique local solution <p with respect to the complex interpolation spaces V = 

We are now going to show that the mappings E, G, Y and Z satisfy Assumption 6.1. 
To that end we need the following preparatory lemma. 

6.6 Lemma. If ^ : R — > R is continuously differentiable, then £ induces a Nemyckii 
operator from L°° into itself which is boundedly Lipschitzian. If £ : R — > R is twice 
continuously differentiable, then it induces a Nemyckii operator from W 1,q into itself 
which is boundedly Lipschitzian. 

The proof is straightforward. Recall that, according to Definition 3.4, q is fixed and 
larger than two. 

6.7 Lemma. Let the Assumptions 3.7, 3.13 and 3.16 be satisfied. Then the equation 
(5.25) defines mappings : [T ,Ti] x L°° — > L°°, k = 1,2, and the restriction of 
each Qk to [T ,Ti] x W 1,q takes its values in W 1,q . Moreover, there is a number 
rj g]0, 1] and then for any bounded subset M C L°° a positive number Qm exists 
such that for all t, t £ [To, Ti] and all if), ip £ M: 

\\Q k (t, ip) " Qk(i, VOHloo < Q M (\t - t|" + \\if} - ^|| L »), k = 1,2. 

Analogously, for each bounded subset M C W 1 ' 9 there is a positive number Qm such 
that for all t, t £ [To, TJ and all if), ip £ M: 

\\Qk(t,ip)-Q k (i,iP)\\ w i, q < Q M (|i-^+||^-^|| wl ,,), k = l,2. 
The proof is obtained from Corollary 5.4. 

6.8 Lemma. Let the Assumptions 3.7, 3.13 and 3.16 be satisfied. If £ : R — > R is 
continuously differentiable, then £ induces operators 

[T , Ti] x L°° 9 (t, ^— > £(Q fc (i, VO) e T°°, fc = 1, 2. 

Moreover, there is a constant r\ £]0, 1] and for any bounded set M C L°° a constant 
£m such that for all t, t £ [To, Ti] and a/i if>, ip E M: 

u(Q k (t, vo) - iu~ < mi* - *p + - $| L ~), ^ = 1,2. 

//£ is twice continuously differentiable, then the restriction of^oQ k to [To,Ti] xW 1 '' 
maps into W 1,q , k = 1,2. Moreover, there is a number rj £]0, 1] and for any bounded 
subset M C W 1 ' 9 a constant £m such that for all t, t £ [T , TJ and a// if}, if) £ M; 

ne(Q*(*, vo) - ik^ < mi* - *p + 11^ - ^ii w m), fc = 1,2. 
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The proof follows from Lemma 6.6 and Lemma 6.7. 

6.9 Lemma. Let the Assumptions 3.7, 3.13 and 3.16 be satisfied. Then there is 
a number rj g]0, 1] such that the mappings E and G defined by (5.28) and (5.34) 
satisfy the conditions (6.1), (6.5), (6.7), and (6.2), (6.6), (6.8), respectively. 

Proof. The functions -kr are continuously different iable by Assumption 2.3. Conse- 



quently, by Lemma 6.8 the mappings E k , given by 

[T j Ti] x L°° 3 (t, if)) i — > — ^4 -GL°°, fc = 1,2, 

y ' F k (Q k (t,iP)) 

are well defined. Moreover, Lemma 6.8 provides a constant t] e]0, 1] such that for 
any bounded set M C L°° a constant Cm exists such that for all t, t e [To, Ti] and 
all if), ^ e M: 

\\E k (t, if)) - E k (i, if)) || L =o < C M (|t - + ||V - ^||l»), k = 1, 2. 

Since W 1,,? embeds continuously into L°° for any bounded set M C W 1 ' 9 there is a 
constant, again named Cm, such that for all t, t G [T ,Ti] and all if), if) E M: 

\\E k {t, if)) - E k (t, ^)|| L =o < C M (|t - £|" + ||^ - ^|| w m), fc = 1,2. 

The identity = — -Efc and Assumption 2.1 now imply (6.1) and (6.7). According 
to Lemma 6.7 the sets 

{Q k (t,4) : (t,<f))e[T ,T 1 }xM}, k = 1, 2, 

are bounded in L°°. Since the derivative of the carrier distribution functions T k , 
k = 1,2, are continuous and positive, (6.5) immediately follows. 

Using the second assertion of Lemma 6.8 we verify (6.2), (6.6), and (6.8) in a similar 
manner. □ 

6.10 Lemma. Let the Assumptions 3.7, 3.13, and 3.16 be satisfied. Then the map- 
ping Z given by (5.27) defines a family {Z(t, i))}( t ^) & [T ,T 1 ]xw 1 ' q of linear, compact 
operators Z(t,<f>) : If — > L p . Additionally, there is a Holder exponent rj e]0, 1] and 
constants Zm such that (6.3) and (6.9) are satisfied. 

Proof. It suffices to show the analogous assertions for the entries of the operator 
matrices Z(t,ifi). Firstly, Lemma 6.8 gives us the estimate 

\\T' k (Q k (t,if))) -r k (Q k (i,i)))\\ B[LP) 

< K(QkM))-rt(QkM))\\L~ 

<CM(\t-i\^+\\iP-iP\\ w i, q ), k = l,2, 
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where the constant Cm can be taken uniformly with respect to t, t £ [7o,7\] and 
■0, ip from any bounded set M C W 1 ' 9 . This estimate together with Assumption 2.1 
implies (6.9). As {Pq 1 ^ is a linear and even compact operator from LP into itself, 
this gives (6.3). □ 

6.11 Lemma. Let the Assumptions 3.6, 3.7, 3.13, and 3.16 be satisfied. Then the 
mapping Y defined by (5.30) meets the conditions (6.4) and (6.10). 

Proof. At first one deduces from the assumptions and Corollary 5.4 that (5.29) 
defines a mapping R : [T ,Ti] x W 1 ' 9 -> L p for which there is a Holder exponent 
r\ e]0, 1]. Moreover, for any bounded set M C W 1 ' 9 exists a constant Cm such that 
for all t, i e [T , 7\] and all ip, $ e M: 

\\R(t,i/>) - R(i,#)\\v < C M {\t - t\ v + 11^ - ^llw^)- 
Applying Lemma 6.9 and Lemma 6.10 one obtains (6.4) and (6.10) for the mapping 
[70,71] x W 1 ' 9 3 (t,if>) .— [J + Z(t,V)]S(*,^)i2(*,V). 

The addends 6^. and i<p' of (5.31) have the required properties due to Assumption 3.7 
and Assumption 3.13, respectively. For V§ x d' they follow from Assumption 3.13 (see 
also Remark 3.15), Assumption 3.16 and the fact that Vo is an isomorphism from 
W~> q onto WT 1 ' 9 . The addend P ^^W¥m o f (5-31) can be treated by means of 
Lemma 6.8 and Assumption 2.1. □ 

We are now going to establish existence and uniqueness of a local solution to the 
evolution equation (5.35). 

6.12 Theorem. Under the Assumptions 3.6, 3.7, 3.8, 3.11, 3.13 and 3.16 the quasi- 
linear parabolic equation (5.35) with the initial condition 4>(Tq) = $° admits a unique 
local solution in the sense of Definition 6.2 with respect to the interpolation space 

V=[h^V]g. 

Proof. According to the Lemmas 6.9, 6.10, 6.11 the mappings E, G, Z, and Y, 
defined by (5.28), (5.34), (5.27), and (5.30), respectively, fulfill Assumption 6.1. 
Hence, the result follows from Proposition 6.5, see also Remarks 6.3 and 6.4. □ 

7 Main result 

We are going to show that a solution of the evolution equation (5.35) in the sense 
of Definition 6.2 provides a solution of the van Roosbroeck system in the sense of 
Definition 4.1. 

We start with a technical lemma. 



Preprint 1189, Weierstrass Institute for Applied Analysis and Stochastics, Berlin 2006 



28 



H.-Chr. Kaiser, H. Neidhardt, J. Rehberg 



7.1 Lemma. Let £ : R — > W be twice continuously differentiable. The composition 
£o-0 is from C([Tq,T], L°°), if i(jeC([Tq,T], L°°) . Ifip composed with the embedding 
L° 0, -^L P ', p > 1, is continuously differentiable in L p on ]Tq,T[, then £ o ip com- 
posed with the same embedding is continuously differentiable in LP on ]To,T[ and its 
derivative is given by 



dt 



■(t) = Z'ty(t))il>'(t)etf, te]T ,T[. 



Proof. If hi, hi G then, by Lemma 5.1 — see also Assumption 2.3, we may 
write 

Z(hi) ~ = £'(/ii)(/n - h 2 ) + T{h x , h 2 )((h 1 - h 2 ) 

where T(hi, h 2 ) converges to zero in L°° if hi^L°° is fixed and h 2 approaches h\ in 
the L°°-norm. Now we set h\ = ip(t) and h 2 = and divide both sides by t — t. In 
the limit i -> t there is \im { _ t T(ip(t), = in while lim^ Ml = ^'(t) 
in LP by supposition. □ 

Our next aim is to justify formula (5.23). 

7.2 Lemma. Let the Assumptions 3.7, 3.8, 3.13, and 3.16 be satisfied and assume 
that <p is a solution of (5.35). We define 

z = (z u z 2 ) with z k (t) = ^ k ( t )+b k (t)+(-l) k i < Po (t), k = 1, 2, t G [T , T], (7.1) 

and (pit) = C(d(t), z(t)) . Then Q k (t, </>(£)) = Zk(t) + (— l) fc jy?(t), and the functions 

[T ,T] 3 t i— > G k (t,$(t)) = G k (QkM(t))) G 

and 

are continuous and concatenated with the embedding L°°^LP they are continuously 
differentiable on ]T ,T[. The time derivative of u k is given by 

u' k (t) = p' k (t)F k (Q k (tJ(t))) 

+ Pk (t)H(Qk(t, ?(*))) [?*(*) + + (-1) Vo(*) + (-1) V(*)] (7.2) 
fc = l,2,te]T ,r]. 

Proof. Due to Assumption 3.8 and Definition 6.2 the function belongs to the space 

C , ([T ,T],L oo )nC ,1 (]T ,T[,L p ) (7.3) 

see also Remark 3.12. Hence, the Assumptions 3.7 and 3.13 ensure that the func- 
tion z also belongs to this space, and by Corollary 5.4, so does the function ip = 
C(d(t),z(t)). Thus, we may apply Lemma 7.1. □ 
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7.3 Remark. Lemma 7.2 justifies the formal manipulations in §5.2. First, (5.23) is 
given a strict sense. Furthermore, the differentiation of Poisson's equation (5.19) has 
the following precise meaning: since <fi is from the space (7.3), the function t i— > ip(t) 
is differentiable — even in a much 'better' space than <p — see Corollary 5.4. Hence, 
the right hand side of (4.4) is differentiable with respect to time in the space ,q 
and (5.19) is an equation in the space W~ ,q . 



We come now to the main results of this paper. 



7.4 Theorem. Under the Assumptions 3.6, 3.7, 3.8, 3.11, 3.13, and 3.16 van 
Roosbroeck's system with initial condition 4>{T ) = <J>° G W 1 ' 9 admits a unique local 
in time solution in the sense of Definition 4-1- 



Proof. By Theorem 6.12 the auxiliary evolution equation (5.35) admits — in the 
sense of Definition 6.2 — a unique local solution <fi satisfying the initial condition 
4>(Tq) = $°. Let us show that — in the sense of Definition 4.1 — the pair {</>,</>}, 
with tp given by 

<p(t) = cp (t)+£(d(t),z(t)), t e [T ,T], (7.4) 

and z according to (7.1), is a local solution of van Roosbroeck's system. First, (4.3) 
is identical with (6.12). By the embedding V Wp' L°° (see Remark 3.12) the 
function [T ,T] 3 1 1— > <p(t) G L°° is continuous, and so is the function [T ,T] 3 t \— > 
$(t) G L°° in view of Assumption 3.8. Thus, G C([T , T], L°°) n C rl QT ,T[,L p ). 
Moreover, for z, see (7.1), one obtains from the Assumptions 3.7 and 3.13 that 
z G C([T , T], L°°) fl C^JTq, T[, L p ). Consequently, property (4.2) follows by Corol- 
lary 5.4, while (4.9) results from Lemma 7.2. The Poisson equation (4.4) with 
densities (4.6) is obviously satisfied by (7.4) due to the definition of C (4.10) fol- 
lows from V0 fc G L7(]T ,T],L 9 ), k = 1,2, and Lemma 7.2. (4.11) is implied by (6.12) 
and (6.13). It remains to show that the continuity equations (4.5) are satisfied. For 
this, one first notes the relations 

Q k (t, 4>{t)) = Mt) + (-l) fc ^(t) + h(t) = z k (t) + (-1) k = 1, 2, (7.5) 

and 

which follows from the definitions (5.25) and (5.29) of R and Q, and (7.1), (7.4). 
Further, in Assumption 3.6 we demand that the mappings r k , k = 1,2, take their 
values in LP — consequently, R takes its values in 17 . From (7.2) and (5.28) one 
gets 
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and by means of the evolution equation (5.35) we obtain 
E(t, 4>(t))u\t) =[I + Z(t, 0(t))] E(t, J(t)) V ■ G(t, ?(t)>V?(t) 

+ [/ + z( t , m] E( t , mR(t, m + ( % d ^m 

We now make use of the representation (4.7) of the currents j = (ji, and get 



E{t,<j>{t)) [u'(t)-V-j(t)-R(t,<f>(t)) 

= Z(t, j(t))E(t, 0(t)) V • j(t) + R(t, cf>(t)) 



We already know that the formal differentiation of Poisson's equation is justified, 
see Remark 7.3. Thus, (5.19) yields 



E{t,<j>{t)) [u'(t)-V-j(t)-R(t,<f>(t)) 

= Z(t, 4>{t))E{t, 0(t)) [V • j(t) + R(t, <P(t)) 

and, observing (5.32) and (7.6), we get 



+ 



u' 2 (t)-X7-j 2 (t)-r 2 (t,v(t),4,(t)) 



0. 



(7.7) 



The operator on the left is continuous on I/; we show now that its kernel is trivial. 
Let /j, / 2 £ L p be such that 



(£) = »■ 



This is equivalent to the relations 



/a 



Ei{t,4>(t)) 
E 2 (t,<t>(t)) 



fi and ^((1 



Ei(tMt)) 

E 2 (t,m) 



-£i(t,0(t))/i 



Vq 1 ^ {[1 + )A) * s a continuous mapping from lV r ' ,? into Indeed, the 

embedding L p ^ W~ 1,q is continuous, and P is an isomorphism between W 7 ^' 9 and 
W* 1 ' 9 , see Proposition 3.3. Hence, we may multiply both sides with fi + pf'~f} fi 



and integrate over Q; this yields 



E 2 (t,cj>(t))- 



Ei(t,4>(t)) 



A fi + 



E 2 (t,<j>(t))' 



fi) dx 



E 1 (t,cj>(t))(l + 



Ei(t,fp)) 
E 2 (t,ct>(t)) 



fidx (7.8) 



Preprint 1189, Weierstrass Institute for Applied Analysis and Stochastics, Berlin 2006 



Classical solutions of drift-diffusion equations 



31 



The quadratic form ip \— > J^(V 1 ip)ip dx is non-negative on L 2 and extends by 
continuity to L p , where it is also non- negative. On the other hand, the function 
Ei(t, 4>(t)) ^1 + J i s almost everywhere on fl strictly positive. Therefore, 

the right hand side of (7.8) can only be non-negative if f\ is zero almost everywhere 
on Q. Hence, (7.7) establishes the continuity equations (4.5). 

To prove uniqueness of a solution of van Roosbroeck's system in the sense of Def- 
inition 4.1 one assures that any solution in the sense of Definition 4.1 procures a 
solution in the sense of Definition 6.2. Indeed this has been done on a formal stage 
by the reformulation of van Roosbroeck's system as a quasi-linear parabolic system 
in §5. In fact, all formal steps can be carried out in the underlying function spaces. 
We accomplish this in the sequel for the crucial points. (4.4) and (4.6) ensure, that 
(f is a solution of (5.4). Hence, Corollary 5.4 implies that </? indeed is continuously 
different iable in W~ q , and, consequently, (5.21) makes sense in W~ 1,q . The deriva- 
tion of (4.6), see also (4.8), is justified by Lemma 7.1. Thus, (5.23) holds in a strict 
sense. The division by Pk^l is allowed because both factors have (uniform) upper 
and lower bounds. The rest of the manipulations up to (5.35) is straight forward to 
justify. □ 



Next we want to establish the natural formulation of the balance laws in van Roos- 
broeck's system in integral form, see (1.1), which is one of the central goals of this 
paper. At first, one realizes that the boundary integral has to be understood in the 
distributional sense — as is well known from Navier-Stokes theory, see [45] — if one 
only knows that the current is a g-summable function and that its divergence is 
p-summable. More precisely, the following proposition holds. 

7.5 Proposition. Let uj C M 2 be any bounded Lipschitz domain. Assume j : uj — > M 2 
to be from L q (uj; M 2 ) and let the divergence (in the sense of distributions) V • j of j 
be p-integrable on u. If q > 2 and p = |, then there is a uniquely determined linear 
continuous functional j„ G W~ l+ ~?' q (du) such that 

j ■ Vip dx + ipV ■ j dx = (j u | ij)\dJ) for allip eW^'iu), (7.9) 

J U) 

1 — — a' 

where (■ \ •) on the right hand side denotes the duality between W •? (du) and 
W~ 1+ 7' q (du) . If, in addition, the function j is continuously differentiate onuo and 
the partial derivatives have continuous extensions to To, then 



j ■ Vip dx + ipV ■ j dx = / ^\quV ■ j da w for all ip £ W 1,q ' {uj), 

J UJ J did 

where v is the outer unit normal of duj, and a u is the arc-measure on duj. 
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Proof. The first statement is a slight generalization, see [30, Lemma 2.4], of well 
known results from [45, Ch. 1]. The second assertion has been proved in [8, Ch. 5.8]. 

□ 

7.6 Theorem. If (<p,(f)) is a solution of van Roosbroeck's system in the sense of 
Definition 4-1, and uj C Q is an open Lipschitz domain, then there are unique 
continuous functions j ku -]To,T] — > W~ 1+ ^' q (du), k = 1,2, such that 

4 / u k (t) dx = (J kv (t) 1 1) + / r k (t, £(<),?(*)) dx, k = 1, 2, (7.10) 



dt . 

where (■ | ■) again denotes the duality between W 1 ~7' q (doj) and W~ l+ ~? ,q (duj) . 

Proof. From (4.5) we obtain for any open Lipschitz domain uj C Q 

u 'k(t) ~ v ' Jk(t) dx = ^ J u k (t) dx - j V ■ j k (t) dx = J r k (t, fi(t), 4>(t)) dx, 

where j k is defined by (4.7). Using Proposition 7.5 we find for every t g]T ,T] 
a unique element j k v(t) <EW~ 1+ ~? ,q (du) such that (7.10) holds. Moreover, continu- 
ity passes over from the functions (4.10) to the mappings ]T ,T] 3 t \— > jkuif) G 
W- l+ 7> q (dcu). □ 

If the currents j k (t) are continuously differentiable on uj and the partial derivatives 
have continuous extensions to uJ, then by the second part of Proposition 7.5 the 
formula (7.10) takes the form (1.1). 



8 Numerics 

Theorem 7.6 is the basis for space discretization of drift-diffusion equations by 
means of the finite volume method (FVM). The FVM was adopted for the numerical 
solution of van Roosbroeck's equations by Gajewski, and this approach has been 
further investigated in [12, 10, 17, 9]. To discretise the spatial domain one uses a 
partition into simplex elements. Let 8 be the set of all edges en = Xi — xi of this 
triangulation, where x±, X2,. ■ ■ are the vertices. Moreover, we define the Voronoi cell 
assigned to a vertex Xi by 

Vi = {x in the spatial simulation domain, such that 

||^ — Xi\\ < \\x — xi\\ for all vertices x\ of the triangulation}, 

where ||-|| refers to the norm in the spatial simulation space M 2 . Now, to get a space 
discrete version of the current-continuity equation, we specify (7.10) with u = Vi, 
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and approximate (jkv(t) | 1) piecewise by jk U a{dVi PI dVi), a being the arc measure 
on the boundary of to = V{. The intermediate value j kil can be obtained as follows: 
The main hypothesis with respect to the discretization of the currents — due to 
Scharfetter and Gummel [49] — is that the electron and hole current density ji and 
ji are constant along simplex edges. This assumption allows to calculate ji a and 
Hn — the constant values on the edge en — in terms of the node values of the 
electrostatic potential and the particle densities, see for instance [17]. Thus, one 
ends up with the following FVM discretization of van Roosbroeck's system for all 
interior Voronoi cells V,: 

e{xi) ^2 (Vy?)izcr(<9Vfc n dV t ) = (d{x t ) + ui(xi) - u 2 (xi)j \Vi\, 
l ■. e it ee 

-Qjr(xi)\Vi\ -j kil a(dVi n dVi) = r k {t, tp, fa, 4> 2 ){xi)\Vi\, 

where \Vi\ is the volume of the Voronoi cells Vi. Here we have tested the Poisson 
equation also with the characteristic function ly. of the Voronoi cell Vi, and we 
have applied Gauss' theorem. In view of Proposition 7.5 we assume, additional to 
Assumption 3.16, d : [To,Xi] — > LP, and observe that ip, can be choosen such that 
(ip, \lvi) = for interior Voronoi cells Vi, see Remark 3.15. Again, we approxi- 
mate the right hand side of (7.9) piecewise by (Vip)ua(dVi fl dV{), and we assume 
- in consonance with the hypothesis about currents — that the gradient of the 
electrostatic potential is constant on the edges of the triangulation, that means 
(V<p)u = {<p{xi) - <p(xi))/\\xi -xi\\. 

Usually, this finite volume discretization of space has been combined with implicit 
time discretization, see for instance [11]. Please note that the strong differentiability 
of the electron and hole density in time is constitutive for this approach. 



9 Outlook to three spatial dimensions 

Much of semiconductor device simulation relies on spatially two-dimensional mod- 
els. However, with increasing complexity of electronic device design spatially three- 
dimensional simulations become ever more important, see for instance [17, 21, 20]. 
This raises the question which of the results for the two-dimensional case carry over 
to the three-dimensional case. In particular, can one expect that in three spatial 
dimensions the divergence of the currents belongs to a Lebesgue space, and is it 
possible to establish strong differentiability of the carrier densities under the rather 
weak assumptions about the reaction terms of this paper. 

Conditio sine qua non for a modus operandi as in this paper is that in the three- 
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dimensional case the operators 

-V • eV : W~' q -> Wf 1,9 and - V • /i fc V : -> Wf 1 '* 

provide isomorphisms for a summability index g > 3. Unfortunately, this is not so for 
arbitrary three-dimensional spatial domains, see [37]. However, one can proof such 
a result for certain classes of three-dimensional material structures and boundary 
conditions, see [7], for instance for layered media and Dirichlet boundary conditions. 
Dauge proved the result in [6] for the Dirichlet Laplacian on a convex polyhedron, 
provided the Dirichlet boundary part is separated from its complement by a finite 
union of line segments. It would be satisfactory to combine this conclusion with a 
heterogeneous material composition. 

Under the hypothesis the afore mentioned isomorphisms exist there are results on 
quasilinear parabolic systems — analogous to Proposition 6.5 — see [43] and [29], 
such that one can obtain classical solutions of the spatially three-dimensional drift- 
diffusion equations very much in the same way as here in the two-dimensional case. 

Acknowledgement. We would like to thank Klaus Gartner for discussions about van 
Roosbroeck's system. 
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